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Abstract 

Using the Bloch representation of an A^-dimensional quantum system and immediate 
results from quantum stochastic calculus, we establish a closed formula for the Bloch 
vector, hence also for the density operator, of a quantum system following a Lindblad 
evolution with selfadjoint Lindblad operators. 
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1 Introduction 



Decoherence is undoubtedly one of tlie main obstacles to use coherent su- 
perpositions of quantum systems as a powerful physical resource, especially 
in the promising quantum processing of information (cf. e.g. [1]). Nowadays 
we may resort to different techniques to deal with this phenomenon [2,3,4], 
each one posing different pros and cons. In this realm, one of the most out- 
standing results is Lindblad-Kossakowski's theorem [5,6], which gives the most 
general form of the evolution equation for the density matrix of an arbitrary 
open quantum system under certain physical assumptions, namely^ trace- 
preservation, complete positivity and Markovianity. The theorem states that 
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the density operator must satisfy an evolution equation with the following 
structure: 



dp{t) 



dt 




(la) 



The operators are known as Lindblad operators; they are completely ar- 
bitrary and H is the (Lamb-shifted [4]) Hamiltonian of the open quantum 
system. Though in the cited works this result was obtained from an axiomatic 
standpoint, other constructive approaches have also arrived at master equa- 
tions with this structure [7,8]. 

In particular, the case when the Lindblad operators are Hermitian plays a 
prominent role in Quantum Measurement Theory [4]. In these circumstances, 
equation (la) can be rewritten as 



This equation has a clear-cut physical meaning: the system on one hand fol- 
lows the unitary dynamics imposed by the Hamiltonian H whereas on the 
other hand it also obeys the non-unitary dynamics represented by the double 
commutators, which tend to project the system onto the eigenspaces of L^. 
When [H, Lk] = 0, the preceding equation represents a typical situation of 
quantum non-demolition measurement. 

On the other hand, another unvaluable tool for studying finite quantum sys- 
tems is their Bloch vector [9,10], i.e. any iV-dimensional quantum system den- 
sity matrix can be expressed as p = ( Jjv + r-X), where Xj (j = 1, . . . , N"^ — 1) 
denote the traceless orthogonal generators^ of SU (N). This Bloch vector r is 
restricted to satisfy certain conditions [9,10]. 

We dedicate this Letter to provide a closed formula for the Bloch vector r{t) 
at an arbitrary time t of an A^-dimensional quantum system satisfying equa- 
tion (lb). Our result rests, apart from the Bloch representation, upon the 
well-known fact that any evolution given by (lb) can be understood as an 
averaged random unitary evolution (cf. e.g. [11]). 

Before presenting the main results, we will clarify the viewpoint followed in this 
Letter in comparison especially with well-known and fully developed stochas- 
tic methods in Hilbert space. There exist two main ways of understanding 
the use of stochatic methods to study the evolution of quantum systems. 
Firstly, the 'quantum trajectory' method (or quantum jump approach) (cf. 
[12,4] and numerous references therein) is profusely employed in Quantum 

^ Those elements satistying (i) A^; = A|., (ii) TrA^ = and (iii) TrA^A^ = 2Skn- See 
[9] and references therein for details. 
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Optics to simulate the evolution of a single open quantum system (normally 
an ion subjected to the action of a laser field) and can be understood as a 
adequate combination of both the standard unitary (Schrodinger) evolution 
and the projective (generalized or not) measurements. The adequacy of this 
technique appears as a consequence of the current experimental feasibility of 
monitorizing the evolution of a single quantum system. 

Secondly, a more fundamental use of these stochastic methods can be seen in 
the several models of dynamical reduction, especially designed to attack the 
measurement problem ([13,14] and references therein). Now the evolution is 
treated as fundamental, with the consequent conceptual deviation from stan- 
dard Quantum Mechanics. 

Here we remain close to the former, but instead of pursuing a physically moti- 
vated simulation of the evolution of a quantum system, we try to find simple 
methods to obtain the analytical solutions of previously known master equa- 
tions. In this sense, the reader may consult [19] where a very simple method 
was presented to solve the phase-damping master equation for any arbitrary 
system. This establishes an important different with the former approach, 
where numerical simulations have been fully developed (cf. e.g. [15]). We thus 
focus only on the mathematical advantages of these stochastic methods, trying 
to identify those unravellings (stochastic evolutions) which allow us to obtain 
the seeked solutions in a simple manner or at least good enough approxima- 
tions to them. 

The paper is organized as follows. In section 2 we briefly revise our needs of 
both the Bloch representation of an A^-dimensional quantum system and of 
stochastic evolution in Hilbert space. In section 3 we use these tools to estab- 
lish our result. Examples are worked out in section 4. Some conclusions are 
included in section 5. 



2 Bloch representation and Stochastic Evolution in Hilbert Space 

2.1 Bloch representation of finite quantum systems 

Given the density matrix p of an arbitrary AT-dimensional quantum system, 
one can always decompose p using the traceless orthogonal generators of the 
Lie group SU {N) and thus write 

P = ^{lN + r-X) (2) 

with r G M^^^"*^. With this decomposition, we can assure that trp = 1 and p"^ = 
p, but some conditions must be imposed upon r to ensure that p > 0, as it must 
be. These conditions were found in [9] and [10] and give rise to the so-called 
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Bloch-representation space ~^), i.e. those [N'^ — l)-dimensional vectors 

valid to represent a quantum density matrix. In mathematical language this 
means that there exists a bijection (f) from the set of density matrix £_|_,i(7Yiv) 
to B{M.^ ^^). This bijective apphcation allows us to study the evolution U{t) 
of a quantum system by focusing on the evolution of its Bloch vector. This 
may be done by forcing the commutativity of the following diagram: 



B{ 



U{t) 



A(t) 



N 



Bi 



) 



i.e. by imposing U{t) — (f)~^ o A(t) o 0. In this sense, once the Bloch vector at 
an arbitrary time r{t) — A(t)[r(0)] is known, we also know the corresponding 
evolved density matrix. 



2.2 Stochastic Evolution in Hilbert Space 



One of the mathematical tools to analyze the evolution of open quantum sys- 
tems in noisy environments is the random evolution in Hilbert space. Though 
there is an impressive mathematical theory behind this approach [16], we will 
only use the following results. An open system obeying equation (lb) can be 
represented by a random state vector which evolves according to a stochatic 
differential equation such as 



d\ct>t) ^ dUM 



(3) 



where Ut is a random evolution operator given by 



Ut = Texp -itH 



k / 



(4) 



where T denotes the time-ordering opeator and B\J(t) denotes the Wa. com- 
ponent of a multidimensional standard Brownian motion [17]. Under these 
conditions, the density operator p{t) at an arbitrary time t > is given by 



p(t)=E[MO)C//] 



(5) 
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In the following we will work out a formula for the Bloch vector r{t) at any 
time t >0 using these results. 



3 A closed formula for the evolved Bloch vector 

Notice that by dropping out the condition of unit trace and positivity, the 
Bloch decomposition can also be applied to any Hermitian operator, say L: 

L^klN + l-X= E ^M^M (6) 

where Aq = In- Note that Iq necessarily satisfies lo = ^ and I G M^^~^ is 
arbitrary (up to the specific details of the operator L) . The generators are 
related through the structure constants of the Lie algebra su{N) [9]: 

[\, Aj] = E ^-ijkh with Cijk = '^■ifijk (7) 

k=i 

where fijk denotes the structure constants of su{N). 

For future purposes we define a multiplication rule between elements of 

given by 

Definition. Let a,P E R'^. The ©-product of a and /?, denoted by a© /3, is a new 
element of M.-'^ whose kth component is given by 

N N 

{a Q 13)^^1 E Cijkaipj (= -2 ^ fijkaiPj) (8) 

i,j=l i,j=l 

It is clear that (M^, ©) has a Lie-algebra structure. We also introduce a second 
multiplication rule by 

Definition. Let Q!,/3 e M^. The H-product of a; and /3, denoted by q;H/3, is defined 

by 

a □ /? = a © (a © /3) (9) 
It is also convenient to define the T-exponential of © and □. 
Definition. Let a, P e R^, with a possibly time-dependent a — a{t). Then we 
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define the deterministic T-exponential of by 



Te"W®/3 = E -T / dh--- dtnTMtn) («(tn-i) (. . . {a{ti) 

(10) 

Equivalent ly for the □ product: 



„=o ^' 

(11) 



To convince oneself that these operations are well-defined, one must notice 
that 



||«©/3||<(C7V)^/^||a||||/3||, (12) 

where C = supj^-^ \cijk\, implies that the series (10) and (11) are convergent. 
These definitions can also be extended to the case in which Ito integrals come 
into play. 

Definition. Let a, P E R^, with a possibly time-dependent a — a{t). Then we 
define the stochatic T^-exponential of © by 



OO 1 ,.f 

T^e"^*)®/? = E -T / dBt, ■■■ fl!B,„T[(a(t,)0(Q;(t,_i)©(. . . (a(ti)0/?)] 

(13) 

Equivalently for the □ product: 



n=0 "^0 

(14) 



Wc will focus on the Heisenberg representation of the density operator p^{t). 
Wc will first give the result for the case of a single Lindblad operator. Equation 
(lb) then reads 

^ = -i[L-(t),[L-(t),p-(t)]] (15) 
where L^{t) — e'*^Le~**^ denotes the Lindblad operator in Heisenberg pic- 



6 



ture. Note that following the preceding decompositions we can write: 



JV2-1 

^=Y1 hi^>^iu = hoXo + h ■ X (16a) 
^i=o 

^= E l^.X^^ = loXo + l■\ (16b) 

/i=0 



For immediate purposes we propose the following 

Definition. Let a e M^^~^ and H — h^X^. Then we define the Heisenberg 

transform of a by 

a^{t) = Te^®a where h = {hi, ■■■ , Hn^-i) (17) 

Note that h is time- independent, so that one can write a^{t) — exp {th(D) a. 
Then we can prove the next 

Theorem. Let p^{t) be the solution of (15) with decomposition p^{t) = {l^ + r^{t) ■ X^. 
Then r^{t) is given by 

r^(t) = E[T,e-^"(*)®r^(0)] = Te^^"(*)°r^(0) (18) 
Proof. Using appendix A we can rewrite (15) as 

p^(t) =E[Te-'/o[^''(^)'-l'^^»][p(0)] (19) 
First apply decompositions (16a) and (16b) to calculate L^it): 

^h + S^-^^-\i-x] = 

= Zo + {Te^®l) ■X = lo + l"{t) ■ X (20) 



f211 



Now we substitute last relation and (2) in (19) to write 

p^(t) = E [Te-<-l^'-]] ^ {In + r^(0) ■ A) 

where for brevity we have denoted = J^l^ {s)dBs. Now the linearity 
of the operator and the stochastic expectaction value yields 

p^(t) = ^ (liv + E[T,e-'"(*)®r^(0)] • A) (22) 

which proves the desired result. The last equality in (18) arises as a 
matter of standard computation of the Gaussian expectation value. 
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A straightforward generalization to nonMarkovian evolutions of this result can 
be immediately obtained: 

Theorem. Let p^{t) be the solution of the nonMarkovian master equation 

^^-:^[L"{t),[L-{t),p-{m (23) 

Let p^{t) = jj (In + r^{t) ■ Xj be its Bloch decomposition. Then r^(t) 
is given by 

r"{t) = E[T,e-'"(*)®r^(0)] = Te^ii"ithimr"(0) (24) 
Proof. Elementary generalization of the preceding proof. 



4 Examples 



We include two simple examples of the use of formula (18). These examples 
are included only to illustrate its possible uses. A detailed exact analytical 
treatment via master equations of these particular examples can be found e.g. 
in [4], chapter 3. For simplicity's sake we will focus on two-level systems, i.e. 
Afc e su(2), with the choice fijk = \^ijk (thus = — x). First we will find 
the Bloch vector of equation (lb) with Hamiltonian H = Uoaz and Lindblad 
operator L = ^/ja^, where cr^ denotes the corresponding Pauli matrix. This 
is an easy example, since both operators commute. The Heisenberg transform 
of I — (0,0,^) is l^{t) — e'^^'^l — as expected. The H-product can be 
represented by a matrix operation 



mr-^(O) = -7Mr^(0) (25) 
where M — (oio). Then r-^(t) can be easily calculated: 

r^{t) = e-^^V^(O) = (e-¥x^(0),e-^y^(0),^^(0)) (26) 

The loss of coherence is clear. As a second example, let us consider H — uqUz 
and L = ^o^- Now they do not commute. The Heisenberg transform of 
/ = (^77,0,0) is calculated through the relation l^ii) — e*^i, where now 

/ u)o \ 

M — I -wo : 

0/ 

^^(^) = V7(cos u)Qt, - sin uj^t, 0) (27) 
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The H-product is implemented through the matrix multiphcation l^{t) □ 
r^(0) = -7iV(t)r^(0), where 



N{t) 



sin^(a'o^) cos cuot sm ivot 
cos u>ot sin LOot cos^ LOot 
1 



(28) 



Taking into account that the integrand T[/^(t„)H(/^(t„_i)H. . . I" {ti)Br" {0))] 
in (18) is completely symmetric, we can write 



This formula allows us to compute r^{t) to any degree of approximation (this 
is especially easy with nowadays computer programs with symbolic compu- 
tation routines). In the present case, the third component can be explicitly 
calculated: z^(t) — e~^z^{0). For the two first components we include the 



computation to second order in 



1 - ^ — ^fl -smc2u;o^) + 



+ 



2U0 2 

-smc uJot + — 



+ 



2uJo 2 
_ 7 (jJpt 
2uJo 2 



7 Wni^ . 2 1 
-smc (Jot + — 



16 V2a;o 



7 \ / 2 2 

— 1 {l + 2ujQt — cos2ujot — 2uJotsm2iOot) 
16 \2tOQj 

2 

{icoot — 2uot COS 2LUQt — sin 2u!Qt 



1 / 7 \ 

^1 -|- sinc2a;o^) H ( ) (1 + 2ujlt'^ — cos 2uot + 2c<;o^ sin 2LOot) 

16 \2u)qJ 



2uq 2 



16 \2ujQ 



{2u)Qt cos 2u!ot — sin 2u!Qt 



x^(0) 



x^(0) + 
(30a) 

(30b) 



where sincx = These relations can be e.g. used to exphcitly check the 
loss of quantum coherence. Working at first order in ^ and choosing r^(0) = 
(1,0,0) we find 



1 - ^(1 -sinc2a;o^) < 1 



(31) 



valid for small times, i.e. 7^ <^ 1. 
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5 Conclusions 



One of the benefits of formula (24) is the possibility of performing a pertur- 
bation approach to the solution of equation (lb). If we explicitly introduce a 
parameter 7 measuring the strength of the Lindblad operator L -w ^L, then 
from (24) one can develop into power series of 7: 



Everything has thus been reduced to a matter of simple computation. Note 
that for the important case of = 2, i.e. of qubits, the abstract operations 
and □ are respectively the vector product x and nested vector product in 'M?. 
Several open questions remain for future work (in progress). Firstly an adap- 
tation of the preceding theorems to infinite-dimensional systems should be 
found. This obviously requires a generalized Bloch decomposition for such 
systems. Secondly an extension of these results to generic Lindblad equations 
is highly desirable. The nonHermiticity of Lindblad operators poses a funda- 
mental obstacle to apply the preceding ideas, since that Hermiticity allows 
us to resort to linear random operators. In the non-Hermitean case, one is 
obliged to use nonlinear stochastic differential equations [18]. 
The stochastic methods have recently been proven to be of great help in com- 
putational tasks regarding some Lindblad equations [19], thus one may right- 
fully hope that a further analysis of these methods will open new paths to 
obtaining the solutions to these master equations. 
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A Computation of E [T exp 



Here we include the computation of the stochastic average of the unitary 

evolution operator given by Texp with a random superoperator 2t(t) = 

/o/5(s)S(s)-rfB„ where /?(s) is a deterministic function, S(s) = ((5i(s),-- - ,©jv(s)) 
is an arbitrary superoperator and denotes multidimensional standard com- 
plex or real Brownian motion. By definition of the T product this is equivalent 
to find the expectation value of 



OO TV It t 

(A.l) 

Using well-known properties of the Brownian motion [17], the expectation 
value reduces to 



n=Oji-jn ^ 



f 

Jo 



dtnT 



(A.2) 
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